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SUMMARY 

Conformal mapping methods are applied to the calculation of 
the effect of channel (two-dimensional tunnel) walls on the ideal 
flow past arbitrary airfoils situated anywhere within the channel. 
The walls of the channel need not be plane but may have any shape. 
The results are compared in specific cases with those obtained by 
two approximate methods, of which the first is a first-order treat- 
ment using image vortices and doublets and the second is a higher- 
order correction developed by Goldstein. 


INTRODUCTION 

In reference 1 a conformal mapping method wa's developed whereby 
the zero- lift velocity distribution could be found for a symmetrical 
airfoil symmetrically located in a plane-walled channel. The purpose 
of the present paper is to extend the previous investigation to the 
case of an arbitrary airfoil situated anywhere within ah arbitrarily 
shaped channel (two-dimensional tunnel). 

The Cartesian mapping function (CMF), Introduced in reference 2 
and used in the method of reference 1, is also used for the problem 
of the present paper. The velocity at any point on the airfoil in 
the channel is found in terms of the CMF and the known conformal 
transformation of a flat plate in a channel. The difference between 
this velocity and the velocity at the same point on the isolated 
airfoil at the same angle of attack represents the effect of the 
channel walls. In order to obtain the velocity distribution on the 
airfoil within the channel, the CMF is applied to doubly connected 
regions analogously to the manner in which Theodorsen's mapping 
function is applied in reference 3. 

The method is given, illustrated numerically by examples, and 
compared with corresponding results by the first-order image theory 
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and by the second-order Image theory of Goldstein (reference 4). 

In addition to the velocity on the airfoil, the velocity on the chan- 
nel vails is obtained by the conformal mapping method. 

SYMBOLS 

The more Important symbols used in this paper are listed as 
follows: 

c chord of airfoil 

Cj section lift coefficient for Isolated airfoil 
Cj ' section lift coefficient for airfoil in channel 
h distance between channel walls 

t thickness of airfoil 

V undisturbed velocity at great distance from airfoil 

v c ^_ velocity on surface of airfoil in channel 

v c2 velocity on chhnnel walls 

v A velocity on Isolated airfoil 

dv velocity correction, v Q ^ - v^ 

a angle of attack of airfoil 

a angle of attack of flat plate 

effective angle of attack of airfoil with respect to curved 
stream 

£ plane of straight lines 

z physical plane 

p circle plane 
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METHOD OF CONFORMAL MAPPING 
The CMF for One Contour 


In previous applications of the conformal mapping method used 
in the present paper (for example, references 1 and 2), a single 
contour such as an airfoil in the physical plane (z-plane) vas 
transformed into a single straight line contour (airfoil chord) in 
another plane (£-plane). The line in the £ -plane is related to a 
circle in a third plane, the p-plane, by a known transformation 
that maps the unit circle vith Its center at the origin into the 
straight line such that the region outside the circle is mapped 
into the region outside the straight line. Because the contour in 
the z-plane also transforms into the same circle in the p-plane in 
such a manner that the regions exterior to the contours correspond, 
the function z - £ is regular everywhere on and outside the circle 
in the p-plane. This vector difference z - £ between conformally 
related points is called the Cartesian mapping function (CMF). 

The real and imaginary parts of the CMF are denoted by Ax 
and Ay, respectively. Because of the regularity of the CMF outside 
the circle, 

CO 

z - £ = Ax (p/p) + iAy (p/p) = C_ n p- n (1) 

0 

where 

p b pe^ 5 


C_n - a- n + ^-n 

On the circle p » e 1 *^ the following relations hold: 

\2li 

Ay (1,<P) = jzl Ax (1,CP*) cot dtp* 


Ax (1,<P) = - -ar/ Ay (1,<P') cot I 


( 2 ) 


Zit 


dtp* 


Equations (2) are the fundamental equations whereby the tr ansf orms - 
tion between the z- and £ -planes can be calculated. 


The CMF for Two Contours 

In general, two contours in the physical plane can be trans- 
formed into two straight-line contours in the £ -plane. The lines 
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in the t- plane can, in turn, be transformed into two concentric 
circles in the p- plane, whose centers are at the origin and whose 
radii are equal to 1 and q (q<l). The transformation is such 
that the region between the two circles is transformed into the 
region between the contours. 

In the case discussed in the present paper one of the contours 
in the z-plane is the airfoil itself; the other contour consists of 
the channel walls, both walls together being considered as one con- 
tour extending to infinity in two directions. The contours in the 
£ -plane consist of a finite straight line into which the airfoil 
is transformed and a transformed channel whose walls are plane and 
parallel to the real axis. In the p- plane, the finite straight 
line, and hence the airfoil, are mapped into the outer circle whose 
radius is unity, and the channels of both the £- and the z-planes 
are mapped into the inner circle whose radius is q. Thus, as the 
outer circle is traced in a counterclockwise direction, the airfoil 
and the finite straight line are traced in a clockwise direction. 

In the same m an ne r, as the inner circle is traced counterclockwise, 
the channel is traced clockwise. 

As in the case of the single contour, the regions at infinity 
in the z- and £ -planes correspond, but the vector difference z - £ 
is regular on the boundary of both circles and within the annulus 
formed by than. As before, z - £ is the CMF and Ax and Ay are 
its real and imaginary parts. Because of the regularity of z - £ 
in the annulus, the CMF may be expanded as follows (cf . equation (l)) 

1 00 

z - t = Ax (p,cp) + iAy (p,cp) C n p n (3) 

-00 

where 


c n - ^ + ib n 

Inasmuch as the full Laurent series is used in equation (3), the 
relations between Ax and Ay on the two circles differ from the 
simple relations given by equations (2). 

Appendix C of reference 2 provides relations between the com- 
ponents of the CMF on the two circles, but the expressions are not 
easily used for the purpose of calculation. More convenient rela- 
tions have been derived in reference 3. Although the correct result 
is obtained, the method of derivation is not fully given. The rela- 
tions are derived in more detail in appendix A of the present paper. 
These relations between the components of the CMF are the following. 
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The subscripts 1 and 2 Indicate the values of the CMF on the 
cirole of unit radius and the circle of radius q, respectively. 
That is: ^ 

Ax x (cp) = Ax(l,<p) 

Ay-^) = Ay(l,q>) 

Aigfa) = Ax(q,cp) 

Ay 2 (cp) * Ay(q,qp) j 


Then, as shown In appendix A: 
AXi(cp) = a 0 + i f A 7l (<P') 


5 cot ^ ^ y 9r> sin n(cp'-cp) 


2 2 
00 


1 1-q 


2n 


dqp'- 


Z f Ay 2 ((p')2^ -^55 sin n(qp'-cp) dcp* 

* J 0 2 1 1-q 211 


1 r 2lt 

AXgCcp) =a 0 - - Ay 2 (cp») 


1 ^ (°R + p sin n(cp'-cp) 

00 


1 1-q' 


dcp' + 


i v” ~n o n ii 

Ay, (<¥>' ) ^ — ^pr Sin n(cp'-cp) dcp' 

* °0 1 1 l_q 2n 


Ay-Jcp) = " n X Ax l^ Cp ' ) || Cot + ^ 2n Sln dcp ' 

; f AXp(cp') ^~^ - 2 % sin n(cp*-cp) dcp’ 

* J 0 1 1-q 

1 

Ay 2 (cp) = b 0 + - AXg (Op* ) 

" f Ax-^qp'jSp ~ 2 \" sin n(cp'-fp) 

* vio 1 1-Q 


i cot IsLfg) + sin n(cp'-cp) 


1 i- c 2n 


dcp* -I 


and 
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The Introduction of elliptic functions simplifies equations (5). 
The elliptic functions introduced at this point and used at other 
places in this paper are treated in various texts with varying nota- 
tions. The notation used throughout this paper is that of Tannery 
and Molk (reference 5). From reference 5 (t. 17, p. 100), the fol- 
lowing series for the ^-functions are obtained: 


1 V*”) X V"' 2n 2n 

* 2 0Ot ''’ + /Z TzH Bl “ (7) 

l . 1-q 


l ^ 4 '(w) V 22 "^ 2 n 

I; sln Zmv < 8 > 


V 
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Hence, from equation (5): 


Axi(cp) 


Ax 2 (cp) 


Ay-^cp) 


Ay 2 (©) 
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In the same way that the relations expressed hy equations (2) 
are a limiting form of Poisson's integral, the integrals in equa- 
tions (5) or (9) are limiting foxms of Villat's analog to Poisson's 
integral. Villat's integral (reference 6) gives the value of a 
function within an annulus vhen the real part of the function is 
known on the bounding circles. 

The relations expressed by equations (5) reduce to those 
expressed by equations (2) vhen the radius of the inner circle 
approaches zero; that is, vhen the channel vails move to infinity. 
The signs differ, however, because the CMF that is defined within 
the annulus in equations (5) is defined within the outer circle as 
the radius of the inner circle goes to zero, whereas, in the case 
of equations (2), the CMF is defined outside that circle. 


The £ -Plane and Its Transformation into the p-Plane 

As already described, the £-plane contains a plane-vailed chan- 
nel within which there is a flat plate. The transformation mapping 
these contours into two concentric circles has been obtained by 
Tamotlka (reference 7), who has also obtained the velocity poten- 
tial for this case. Tamotika's results will be briefly presented 
and the form in which they are most useful in applying the CMF 
method will be given in more detail. 

Let a be the angle of attack of the flat plate. The trans- 
formation between the I-plane and the p-plane is shown in figure 1 
and given mathematically as 


Zh e i(|-fS)y^ aln nCp 2 x 

* ' n(l + 2q 2n cos 28 + q^ n ) 

[pPfe 1 ® + q 2n e _i ®) + p" n (e“ i8 + q 2n e -5 )] +t 


or in another form 


. _ 12h 

^ = it 


A. log + 

21 AOg e _icp + 

p-qe « - 


q 2n sin n<?t> 


j n(l-2q 2n cos 2a + q 4n ) 


4m 



(e 2i “ 



> +T 


( 10 ) 


( 11 ) 


J 
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where h is the distance between the channel walls, 8 = 5 - a, 

icp “ 

and T is a constant. The substitution of p = e in equa- 
tion (10) yields the equation of the flat plate; the substitution 

of p = qe^ in equation (11) yields the equation of the channel 
walls, which are parallel to the real axis. The use of the two 
forms of the transformation simplifies the resulting equations in 
cp for the flat plate and channel walls. 


Four values of the central angle cp (cp^q^opj, ,and<j^) are 
Important in the mapping. From reference 7 the points p = e 
e i(2fl-cpj.)^ denoted by B and B' map into the stagnation points 

. the f lal 
i(2fl-cpg) 


on the flat plate for zero circulation; the points p = qe^ 2 . 


qe “ , denoted by H and H' map into », respectively; the 

icp, m 

points p = e , e , denoted by A' and A, map into the 

extremities of the plate. The points are shown in figure 1. 


The values of q and cp 2 in equations (9) and (10) are deter- 
mined by the length of the plate, its position, and the various rela- 
tions between the four special values of cp. From equation (10) or 
equation (11) 


or 


CP 3 + <P 4 = 2a 
^4 \ 2rt / * 4 \ 2n 


cp 1+ c^ 

*4 VST" 




4\ 2 it 


= 0 


8it 


a n 

^ u 2n cos ncj^ sin ncf^> = 0 


fD 3 _fD 2\ „ / cp 4" cp 2\ . f Cp 3'* tp 2\ . (^4^2 


4\ 2n J v 4\ 2it 


v 4\ 2 it 4\ 2 it 


( 12 ) 

(13) 


(14) 



n n(V z -%) 

q n sin ncf^ sin = 


n(l-2q 2n cos 2a + q 4n ) 


[cos(n-l) a - q^ n cos(n+l) a] 


( 15 ) 
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n n(<P 3 -q> 4 ) 

q sin njpg cos g 


n(l-2q 2n cos 2a + q 4n ) 


[sin(n-l) a - q 2n sin(n+l) 2 ] (16) 

vhere L is tbe length of the plate and d is the vertical distance 
of the midpoint from the center line of the hhannel. 


In order to find the values of q,cp 1 ,cp 2 ,cp 3 , and cp 4 for a 

given length position of the plate, equations (12) to (16) should 
he solved simultaneously. 


In principle, it is possible to transform the flat plate at 
any value of a into the airfoil at angle of attack set at a. 

The value of a is fixed at the value that accomplishes the trans- 
formation with the least labor* In the case to be calculated, <Z 
is set equal to zero. For this value of a, Tamot ilea's formulas 
(reference 7) are considerably simplified. If, the distance between 
the channel walls is taken as unity, the equations simplify as 
follows: 


i 


oo n 

E q sin rfpg 
A ntl-l 2 ”) 


(p“-P-“) . T 


(17) 


or equivalently 




, *?2 
— log e -rz- + 


21 


-iq>„ 


p-qe 



q 211 s in nflpg 
n(l-q 2n ) 






+T (18) 


and «P 2 


can be found by using the equation 


d = 



Jt 


(19) 


The quantities q, <P 3 , and may be found by solving 

simultaneously equations (12), (13), (14), and (15), which also 
become simpler than the equations for the general case. The 
constant t may be restricted to real values, because it merely 
determines the position of the channel and the flat plate with 
respect to the axes in the £ -plane. 
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A special case useful for numerical work is that for which 
d s 0. For this position of the flat plate op 1 = cp g = cpg = -cp 4 = | 
and 


L = 



,2n-l 




( 20 ) 


Hence, q r^n be found from equation (20) alone when the length of 
the plate is prescribed. 

If £ is separated into its real and imaginary parts, 

£ s | + It), | and t) be found as functions of <p. The equa- 
tion of the flat plate is found hy setting p = e iCp in equa- 
tion (17). Then, when the subscript s 1 and 2 denote the values of 
the function on the plate and on the channel, respectively, 


£ 


1 " 



q 11 sin ncpg 

,, 2n* 

n(l-q ) 


sin ncp + t 


( 21 ) 




= 0 


( 22 ) 


Thus the flat plate lies on the real axis of the £-plane. The 

equation of the channel walls is found by setting p = qe iCp in 
equation (18). Then 


t-2 



sin \ (cp-<^) 
sin | (cp4<^) 



q^ n sin nff^ 

nd-q 211 ) 




sin ncp>+ t 




(23) 


^ 1 
^2 = * “2 


d 



(2rt-q> 2 > cp > cp 2 ) 


l (24) 


(CP 2 > cp > - cp 2 ) 
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Airfoil Position and Adjustments in Terms of the CMF 


The z-plane and the £-plane are shown superimposed in figure 2 
in which the geometric meaning of the CMF is also indicated. If 
the abscissas and ordinates of the airfoil are denoted by x^, y^ 

and the abscissas and ordinates of the channel walls by x 2 , 7z> 

the definition of the CMF shows that 


Xj/cp) = li(«P) + Ax x (cp) I 
y 1 (cp) = Ay^cp) 

> 

x 2 ( c P) ■ l 2 (cp) + AXgfa) 
y 2 (fp) = t) 2 + Ay 2 (cp) 


(25) 


In order to determine the constants q and cp 2 that appear 
explicitly in the expressions for | and q and also the angles qpjj 
and ep T that correspond to the leading and trailing edges of the 
airfoil, the airfoil is placed in a normal position with respect to 
the y-axis. If c is the chord of the airfoil and a is the angle 
of attack, the normal position is given by 


= - | cos a 

X 1 ) = I cos a 


(26) 


From equations (25) and (26), the following formula is obtained: 


f x (cp T ) - li(cf^j) = c cos a - Ax^cpij) + Ax^Cffo) (27) 

The angles cp^ and cp T corresponding to leading and trailing edges 
are obtained from the condition of a maximum for the abscissa x^(cp), 


dx^cpjj) 

dcp 


= 0 


dxitoj,) 


or, by equations (20) and (24), 
dAX x (cp) 

n i 

1 


dcp 


_ ? 8 
" *hr f l - 


dcp 


sin ncfg 


= 0 


2n 


cos ncp 


for cp = cp 


N 


or cp T . 


(28) 


(29) 
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The value of cpg, or what is equivalent the value of d, 
found from 

\2n r>2n 


X CJl 

Ay^cp) dcp = Ay 2 (cp) dcp 


Is 

( 6 ) 


as follows: Let r(cp) denote the value of the ordinate of the air- 

foil measured from the center line of the channel . in the t -plane. 
From the definition 

r(cp) = Ay^(cp) + d (30) 

Hence, using equation (6), 

i a 2ji i ~2jt 

f r(cp) dcp - ~ P Ay 2 (cp) dcp (31) 

4) 'V 

and cp2 is obtained from equation (19). 


The constant t is obtained by adding the equations of (26). 
The resulting formula is 


T 



q 11 sin n<?v> 


n (iJi N -rtL) nfPu-cUr) 
sin 0 cos s 


AXi(cpjj) + (cprp ) 

2 

(32) 


These equations completely determine the constants q, T, Cp K , 
cpj, and cp 2 in terms of the CMF. The value of c^ is calculated 
from equations (31) and (19); the use of this value in equations (27) 
and (29) permit these equations to be solved simultaneously for q, 
cpjj, cp T ; and finally T can be determined from equation (32). 


Velocity Distribution on the Airfoil 
and on the Channel Walls 


The complex velocity potential W, derived from the results of 
reference 7, is 


A (l l°8e P+ Cp 2\ ( ^2 \ 
u Vh . 4 \ Zn ) 4 V 2 it ) _ IT 

W = All TZZ ZT7n\ /m \ “ Pir 


/llog e P^ 2 \ 


2« 


2 it l0g e P 


(33) 


2k 


where V is the velocity at infinity and T is the circulation. 



14 


NACA TN 1899 


The velocity distribution on the airfoil and on the channel 
vails is obtained from the velocity potential given by equation (33). 
The formula for the velocity in the z-plane is 


v z s 


d¥ 

dE_ 


dz 

dp 


(34) 


On the airfoil, from equations (8) and (33) 

'/ -S2- „n 


(i 


ie 


dw\ _ 

dp 4=e»“ 2 ” 


-icp 


[ XT' q sin \ „ 

1 8Y ) cos nsp) - T 

v 4 - i-o 2 * ) 


(35) 


\ 1 

The circulation T in equation (35) is adjusted to satisfy the 
Kutta condition at the trailing edge of the airfoil (~) = 0. 


The result is 


00 _ 

E q u sin nq>o 
-E 

* H 


COS 




(36) 


Also 


dz 

dp 


iqp 

p=e ^ 


. (£ \ ,-W . -ltp 

~W) & + dq> e 


41e“ i<p \ ' « aln _ idd^lW _ lm 

' ~~ " 2n 008 “ ,p dip— e * + 

1 1-4 


> (37) 


dAy 1 (cp) 


dcp 


i«P 


J 


Hence, the velocity distribution on the airfoil is 


Y cl _ 4 
V " n 


00 

E 


q 11 sin nfp 0 


i-q 


2 ^ - - (cos n<p- cos irpip) 


( 38 ) 


4 V" q Sin “^2 


n 


“ 2n“ 008 “P - “d^ 

1-q 


dAx x (cp) 

2 

dAy x (qp) 

dcp 

+ 

dcp 


vhere v ^ has been vrltten for v 2 . 



The velocity distribution on the channel walls is found by replacing p by qe x ^ in equa- 
tion (33). The substitution results in 
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Lift on the Airfoil in the Channel 

The lift on the airfoil in the channel can be found by evalu- 
ating a modified form of Blasius' integral in the p-plane. The 
expression for the lift involves the CMF and the radius of the inner 
circle; that is, the lift depends on the shape and the position of 
the airfoil and on the shape of the channel walls as well as on the 
circulation. This dependence is in contrast to the case of the 
isolated airfoil, in which the lift on any body is the same for a 
fixed circulation. The dependence of lift upon the airfoil shape 
for the case of the airfoil in a plane-walled channel has also been 
shown by Havelock (reference 8) who finds the potential function 
directly without the use of conformal mapping. 

The expression for the lift is too complicated for numerical 
calculation. A more convenient way of obtaining the lift is to 
integrate the pressure distribution on the airfoil or the pressure 
distribution on the walls. 


Method of Successive Approximations for Obtaining CMF 

The CMF can now be calculated for a given configuration by a 
method of successive approximation analogous to that of reference 2. 

1. The airfoil and the channel walls Eire drawn such that the 
airfoil is in the normal position, as shown in figure 2. The center 
line of the channel in the £-plane is located on the figure in order 
that the airfoil ordinates r(cp) may be read. The scale is so 
chosen that the distance between the channel walls in the £-plane 

is unity. 

2. From a previous approximation, approximate values of q, 
T,cp 2 /P N , and are known, as well as approximate values of the 
abscissas x^(cp) and x£(cp) at a convenient set of values of cp 
from 0 to 2n radians. Through the use of the known values of x^cp), 
r(cp) is measured. A set of values of Aygfa) are also measured 
through the use of the known values of x 2 (cp). A value of d and 

a new value of cp g are obtained from equations (31) and (19). 

If no better values are available, the initial approximation 
for x x (cp) and x 2 (cp) may be that obtained for the flat plate 
situated along the center line of a p l a n e-walled channel. In this 
case x^ and x 2 are given by equations (21) and (23) for 
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and ® ie value of q is obtained from equation (20), where L 

is replaced by c cos a. Both cpg and 0Pjj equal k/ 2 and CP T 
equals 3«/2. The constant T equals zero. 

3. The functions Ax 1 and Ax 2 are calculated by means of 
the first and second equations of (5). The value of q used is 
the approximate value of step 2. The numerical details of the 
calculation are given in appendix B. 

4. New values of cp N , cp^,, and q are obtained by solving 
equations (27) and (29) simultaneously for these quantities. 

An alternative method of determining cpjj, cp and q is a 
purely graphical one. The approximate function x^(cp), which is 
also a function of .q, is plotted against cp in the regions of 
the extreme values of x^. From this graph cp N and cpj are 
determined. These values are substituted in equation (27), from 
which a new value of q is obtained that is used to re-evaluate x^. 
The procedure is continued until sufficient accuracy is obtained. 
Finally t is calculated from equation (32) . 

5. A new set of values for x^(cp) and XgC^p) are calculated 
using the new values of the constants and the values of Ax^ and Ax 2 
calculated in step 3. 

Steps 2 through 5 are repeated until a plot of y(cp) against 
x(cp) for both the airfoil and the channel walls yield shapes, that 
are as close as desired to the shapes plotted in step 1. 

If the walls of the channel in the z-plane are flat, Ay 2 (cp) 
is set equal to zero, and a considerable simplification in the 
numerical procedure results. This case is the most common and the 
method is not at all difficult to apply. The discussion of numerical 
results will provide an idea of the actual work involved. 

After the components of the CMF and the various constants have 
been evaluated by the method of iteration Just described, the veloc- 
ity distribution may be found from equations (38) and (41) for the 
airfoil and for the channel walls, respectively. The derivatives 
of the CMF in the formulas for the velocity distribution were 
measured in the cases calculated; although an expression exists 
that gives the values of the derivative in terms of the CMF as in 
reference 1, it is too cumbersome to use. 
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ILLUSTRATIVE EXAMPLES USING CONFORMAL MAPPING 

The method of conformal mapping outlined has been applied to 
the 12-percent symmetrical airfoil treated in reference 1. The 
ordinates of this airfoil are given in table I and the airfoil shape 
is shown in the figures in which the velocity distributions are 
plotted. For the calculations of the present paper the airfoil was 
assumed to be placed at the center of a plane-walled channel. The 
chord to height (c^x) ratio was taken to be 0.5. Velocity correc- 
tions were calculated for angles of attack of 0° and 4°. 

For the case of a = 0° the range of <p from 0 to 2 n radians 
was divided into 24 equal intervals. Two approximations, starting 
from the x(op) of the flat plate, were necessary for the derived 
airfoil contour to coincide with the given contour for a scale of 
chord length of 20 inches and ordinate scale five times that of the 
abscissa scale. In no case were more than six terms used in any of 
the infinite series in the preceding formulas, for the series con- 
verge rapidly. The velocity distribution for the case of a = 0° 
is shown in figure 3. The velocity distribution on the walls of the 
channel is included in the figure and is drawn to a scale five times 
as large as the scale for the velocity distribution on the airfoil. 
The CMF together with the velocity distribution is given in table II 
The velocity distribution on the airfoil for this case had been 
previously calc ula ted by the method of finite chord in reference 1. 
The results are compared in figure 4 and are in close agreement, 
which indicates that the numerical methods used in both processes 
were sufficiently accurate. 

The velocity distribution for the case of angle of attack of 4° 
is plotted in figure 5. Figure 6 shows for the purpose of compari- 
son the velocity distribution for the airfoil in the free stream at 
a = 4°. In this case four approximations, starting from the flat 
plate, were necessary to obtain coincidence between the derived air- 
foil and the given airfoil to the same ordinate and abscissa scale 
as in the case of a = 0°. In the first three approximations the 
cp range was divided into 24 equal intervals, but in the fourth 
approximation the length of the intervals was halved so that the 
CMF was evaluated at 48 points. The mapping data and velocity dis- 
tribution are given in table III; the nature of the CMF is shown by 
figure 7 where the component functions are plotted. The velocity 
distribution for the airfoil in the free stream was obtained by the 
method of reference 2. 

The velocity correction for the airfoil at an angle of attack 
of 0° was discussed in reference 1. The velocity corrections for 
the airfoil at the angle of attack of 4° are plotted in figure 8. 
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The irregularities of the correction are due to local curvature 
fluctuations of the airfoil surface and correspond to the irregu- 
larities found in the corrections for the same airfoil at a = 0°. 
(See reference 1.) 

The velocity corrections are positive on the upper surface of 
the airfoil hut are for the most part very nearly zero on the lower 
surface. This behavior of the correction indicates that the lift 
on the airfoil in the channel is greater than that on the airfoil 
in the free stream. The increase in lift has been shown by other 
authors through the use of approximate methods (see references 4, 

7, 8, and 9) and will be further discussed. 

The influence of the airfoil on the velocity distribution on 
the channel walls is shown in figures 3 and 5. The velocity dis- 
tribution on the walls is very sensitive to the angle of attack. 

When the angle of attack is 0° (fig. 3) the nondimens ional velocity 
on both the walls is greater than unity. The velocity rapidly 
approaches unity both upstream and downstream of the airfoil until 
at 1.75 chord lengths upstream and downstream of the origin the 
velocity has decreased from its maximum value 1.03 to substantially 
the value 1. 

In contrast, when the angle of attack is 4° (fig. 5), the 
velocity is less than unity on the lower wall, and on the upper wall 
the velocity markedly increases over the velocity for the case of 
a = 0°. The m»Yimum velocity on the upper wall moves forward toward 
the position at which the airfoil approaches closest to the wall; at 
the same time the m inim um value on the lower wall is located at the 
position near the leading edge where the zero streamline rises to 
meet the airfoil at the stagnation point. On both the upper and 
lower walls the velocity approaches unity less rapidly than in the 
case of a = 0°. On the upper wall the maximum velocity is 1.095; 
the velocity 1.75 chord lengths upstream of the origin is 1.013; the 
velocity 1.75 chord lengths downstream is 1.010. On the lower wall 
the minimum velocity is 0.965; the velocity 1.75 chord lengths both 
upstream and downstream is 0.990. 


APPROXIMATE VELOCITY CORRECTIONS FOR AN AIRFOIL PLACED 
ALONG CENTER LINE OF A PLANE-WALLED CHANNEL 

If an airfoil is placed midway between the walls of a plane- 
walled channel, simple approximate velocity corrections may be 
derived under the conditions that the angle of attack is small and 
that the thickness, chord, and camber are small in comparison with 
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the dimensions of the channel. Two such corrections will he 
explained. Both corrections depend upon the successive reflection 
of the airfoil in the channel walls by which a cascade of airfoils 
alternately upright and inverted is obtained. As is well known 
(see reference 9), the flow through such a cascade is equivalent 
to the flow about the airfoil in the plane-walled channel. In the 
first-order approximate theory, the image airfoils are replaced by 
doublets and by vortices; in the more elaborate treatment developed 
by Goldstein (reference 4), higher-order singularities are included. 
Inasmuch as the method of conformal mapping developed in the present 
paper is applied numerically to a symmetrical airfoil at the center 
of the channel, the approximate theories will be quantitatively 
discussed only for such airfoils. A more general treatment would 
follow along similar lines. 


First-Order Theory 


In the development of the first-order theory the vortex and 
the doublet are assumed to contribute independently to the velocity 
correction. The effect of the image vortices is to curve the stream 
and to increase the effective angle of attack and lift on the air- 
foil in the channel. The image doublets increase the velocity at 
the center of the channel and thus take into account the constricting 
effect of the channel walls. Glauert (reference 9, p. 49) obtained 
a formula for the ratio of the lift in the free stream to the lift 
in the channel. If it is assumed that the vortices merely change 
the angle of attack, the Kutta condition combined with Glautert's 
formula yields the following result: 


sin a it 2 / c V 

sin 24 \ h) 


(42) 


where a is, as before, the angle of attack with respect to the 
direction of the flow at infinity and is the effective angle 
of attack due to the curved stream. 


The increase of velocity at the center of the channel induced 
by the image doublets is assumed to be that due to the airfoil at 
its angle of zero lift. If this increase is denoted by u and, 
as before, V is the velocity at infinity in the channel, the fol- 
lowing relation is true: 



( 43 ) 
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where for symmetrical airfoils 

* - $ Ijr f * (!) <«> 

as in reference 9 (p. 55). Here v c ' is the velocity on the air- 
foil when the airfoil is in the channel at an angle of attack of 0° 
And y ia the distance to the upper surface of the airfoil measured ' 
normally from the chord line. The integral in equation (44) is 
taken with respect to the surface distance s along the upper sur- 
face of the airfoil from leading to trailing edge. 


In the calculation of the strength of a doublet that is to 
replace an isolated airfoil, Vj^ rather than v c ' should be used. 
However, Inasm uch as the strength of the doublet must be increased 
when it is used to replace the same airfoil in cascade, the use of 
v c ', which is greater than v^, will change the value of X in 
the right direction. 


The velocity correction is defined as 


Ay /M (Ji\ 

v 58 w ai ~ W« 


(45) 


where ( -%• ) is the velocity on the airfoil in the channel expressed 
a l 

as a fraction of the ultimate upstream velocity when the airfoil is 


at an effective angle of attack and where is the isolated 

airfoil velocity for the angle of attack a. Since the airfoil is 
small compared with the breadth of the channel, the flow about the 
airfoil in the channel is equivalent to the flow about an airfoil at 
an angle of attack in a free stream whose velocity at a great 

distance away is V + u. Therefore the following relation is true 



The result is that 
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The formula for the velocity correction shows the Importance of the 
changed angle of attack, for one part of the correction is the 
difference in the isolated airfoil velocity distributions at angles 
of attack a and the other term of the correction is propor- 

tional to the isolated velocity distribution at the increased angle 
of attack. 


The correction obtained by the use of vortices and doublets 
is valid to the first order in and tc^i 2 . When the angle 

2 


of attack is 0°, the parameter does not appear (reference 1) 


Goldstein's Second-Order Velocity Correction 

Goldstein (reference 4) first replaces the image airfoils by 
the doublet, the vortex, and the higher-order singularities given by 
the potential function of the airfoil in a uniform free stream. The 
nonuniform disturbance velocity produced by these singularities in 
the physical region, in particular at the location of the physical 
airfoil, is calculated, taking into account the change in direction 
of the stream. This first-approximation nonuniform disturbance 
velocity (a) changes the velocity distribution on the airfoil from 
its isolated free- stream value and (b) changes the value of the 
singularities that are to be imaged. Change (b) is evaluated and 
a second-approximation nonuniform distribution of the airfoil in the 
final nonuniform stream is calculated. 

In principle, Goldstein's method is capable of yielding to any 
degree of accuracy the effect of a plane-walled channel on the two- 
dimensional velocity distribution of an arbitrary airfoil, arbitrar- 
ily situated. The successive approximations become increasingly 
laborious, however, and only the second-approximation formulas are 
given in reference 4. 

The second-approximation formula for the constriction correc- 
tion for the symmetrical airfoil situated in the center of the chan- 
nel at a small angle of attack is obtained as: 

T c / u\ " £(*) + sin (© + ox) + sin a-J 

Vj \V/ [sin (0 + a) + sin a] 


( 49 ) 
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so that 


Av /Jc 

V - V 


(50) 


yhere U here represents the sum of the ultimate upstream velocity 
and the velocity at the center of the channel induced by the singu- 
larities so that ? - 1 corresponds to u/v of the first-order 

theory; is, as in the previous approximate theory, an effective 

angle of attack with respect to the direction of the stream; the 
function p(e) is a measure of the distortion of the stream caused 
by the singularities. 


The Goldstein second-order image correction is accurate to the 


(C\ c tc /I 

A /°\ 

\h)’ h 2> U 

0 ’ U) 


c 3 t 


c 2 t 2 
h 4 J 


ct' 


and 


(*)• 


When the 


angle of attack is zero, the terns 
not appear. 




Discussion of Numerical Besults 
of Approximate Theories 

The first-order and second-order corrections were calculated 
for the 12-percent- thick symmetrical airfoil. The corrections for 
the airfoil at zero lift have been discussed in reference 1. The 
results for the angle of attack of 4° are plotted in figure 8. The 
constants used in the first-order correction are 


X 

u/V 

®1 

3.93 

0.0116 

4.459° 


Those for the second-order correction are 


C 0 

C 1 

C 2 

C 3 

C 4 

“1 

?- i 
V 

0.08722^ 

0.05534 

-0.02401 

0.00455 

0.00475 

4.200° 

0.0108 


The first-order theory yields good results for the upper sur- 
face of the airfoil in that the correction so derived shows the same 
over-all trend as the correction obtained by conformal mapping. The 
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approximate correction appears to be a mean curve to which are added 
components due to the curvature of the airfoil. For the lower sur- 
face, the approximate correction is not quite so good a mean line as 
it is for the upper surface. For both upper and lower surfaces, the 
contribution to the velocity correction due to the doublets and that 
due to the change in angle of attack are equally effective in form- 
ing the total correction. 

For the upper surface of the airfoil, the Goldstein second- 
order image correction follows the same trend as the first-order 
image correction, but the values are more nearly constant. The 
second-order correction for the lower surface follows more closely 
the trend of the mapping correction than the first-order correction. 
From this example, the second-order correction appears to be more 
accurate than the first-order correction. 

The incremental velocities u/v and 2 _ q of the first- M 
second-order corrections, respectively, are in good agreement but 
the values of the effective angles of attack Oq_ differ markedly. 
This difference accounts for the difference in the nature of the 
correction curve of figure 8 near the leading edge of the airfoil. 


CALCULATION OF LIFT AND MOMENT 

For the case of angle of attack of 4°, the lift coefficient c, ' 
for the airfoil in the channel was calculated by integrating the 4 
pressure distribution about the airfoil. The calculation for c^' 
was also carried out by means of the two approximate theories. 

The isolated airfoil lift coefficient c^ was 0.478. The 
value of Cj' obtained by the integration of the pressure distri- 
bution is 0.537; that value obtained from the second-order theory, 
0.522 by the formulas of reference 4; and that value obtained from 
reference 9 (p. 49), 0.532. All the values of c^' obtained indi- 
cate the expected increase in lift for the airfoil in the channel 
and also show good agreement among themselves in that they do not 
vary more than 3 percent. The lift-coefficient correction, c^ ' - cj, 
varies, however, about 30 percent among the different theories. 

The lift coefficient Cj ' was also calculated by integrating 
the pressure distribution on the walls of the channel. Theoretically, 
the integration should be carried out to infinity on either side of 
the airfoil. The practical calculation is, of course, impossible. 
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The integration is therefore carried out only over a finite range 
to yield the lift coefficient Cj", and a correction factor used 
to take into account the effect of the rest of the channel. 

The correction factor tj, which is equal to c has 
been derived in an approximate form in the appendix of reference 10. 
The airfoil is replaced by a row of vortices, which are imaged in 
the walls of the channel. The q factor for an individual vortex is 
calculated. The final q factor is obtained by averaging rj for 
each vortex with a loading derived from thin airfoil theory as a 
weighting factor. 

In figure 9 the lift coefficient c j” is plotted as a function 
of the limits of integration, which were taken symmetrically about 
the origin. The value of cj”, obtained by integrating the pres- 
sure distribution 1.75 chord lengths upstream and downstream of the 
origin, is 0.493. When this value is divided by the value of c^', 
derived by integrating the pressure distribution on the airfoil, a 
value q of 0.918 is obtained. The value of t} obtained by the 
method of reference 10 is 0.900. The value erf c j', obtained from 
the appro xima te value of rj, is 0.548. The correction factor 
obtained by the approximate method is satisfactory to the order of 
the approx ima te theories previously discussed. 

It is also possible to obtain the moment on the airfoil about 
any point by integrating the moment of the pressure (accurately 
calculated) on each element of area on the channel walls. A factor 
analogous to the rj factor can be so determined that the integration 
for the moment over a finite range may be extended to take into 
account the regions on the channel walls a great distance away. 


CONCLUSIONS 

Tha analysis and numerical calculations of the present paper 
lead to the following conclusions: 

1. The method of conformal transformation by means of the 
Cartesian mapping function provides a satisfactory numerical solu- 
tion to the problem of obtaining the local velocity corrections for 
an arbitrary airfoil in a channel for the case of two-dimensional 
frictionless incompressible flow. 

2. If closeness to the velocity corrections obtained by con- 
formal mapping is used as a criterion, the second order Cioldstein 
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correction is more accurate than the first-order image vortex and 
doublet correction for thin airfoils at small angles of attack in 
giving velocity corrections in the examples calculated. 

3. If it is necessary to obtain a higher-order correction than 
the second, the method of the Cartesian mapping function is probably 
more convenient to use than the Goldstein type correction. 

4. The channel lift coefficients obtained by the two approxi- 
mate theories are in good agreement with the lift obtained from the 
mapping velocity distribution; the lift corrections obtained by the 
two approximate theories are not in good agreement with the correc- 
tion obtained by mapping results. 

5. The existing method of finding the lift coefficients from 
the velocity distribution on the channel walls has been satisfac- 
torily checked. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, December 4, 1946. 
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APPENDIX A 


DERIVATION OF THE RELATIONS BETWEEN THE REAL 

AND IMAGINARY PARTS OF THE CMF 

Inasmuch as the CMF z-£ is regular within the annulus and 
also on the hounding circles in the p- plane, it may he expanded in 
a Laurent series, which is valid in the annulus and on the circles 
hounding the annulus. Thus 

03 

z - £ ■ Ax(p,q>) + iAy(p,<?>) = Cnp 11 (Al) 

-00 

From equation (Al) the following expressions are obtained: 

AXj/cp) = a Q + a_ n ) cos ncp - ^(b n - h_ n ) sin nqp 

1 1 


CD 

Ay^fcp) = h Q +22(a n -a_ n ) sin nqp + 2Jb n +b_ n ) cos ncp 
00 00 ^ 

AXg(cp) - a Q +X/(a n q n + a_ n q” n ) cos i£p - 2_J(h n q n - h_ n q _n ) sin nCp 


00 


HA2) 


Ay 2 (cp) = b +2Z(a n ^ n - a. n q' n ) sin nq> + ]T)(b n q n + b_ n q" n ) cos ncp 


The values of a n and h n can be found by means of Fourier's rule 
in terms of the CMF. 

When Sq and b^ are evaluated, the conditions of consistency 

that are necessary conditions for the regularity of the CMF in the 
annulus appear as 


2jta Q = P Ax,(cp) dcpef 2 * AXg(cp) dcp 

p2jt pZit 

2nb 0 =J Ay x (cp) dqp =J Ay 2 (cp) dcp 


(A3) 
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Four relations are desired: Ax^ and Axg expressed in terms 

of Ay-j^ and Ay g and, conversely, Ay^ and Ay g expressed in 
terms of Ax^ and AXg. The derivation of the expression for Ay^ 
in terms of Ax^ and Axg will now he carried out. The other 
relations will follow analogously. 

Through the use of the first and third equations of (A2) 
through the use of Fourier's rule, the coefficients a n and b n are 

evaluated. As a result of the calculation, the following equations 
are obtained: 



The values of the coefficients a Q and b n are substituted in the 
infinite series expression for Ay-j^ given by equations (A2). The 

values for Ko, Kg, ®i# and as Siv® 1 * by equations (A5) are 
also used. Thus 

00 p2it 

1 V — \ n -n 

dy^i 0 ?) =bQ +- ^ ^ n + *^_ n Ax^(<P' )(sin itfPV cos n<p- cos nSP' sin i*p)d<p' + 

1 q _<1 Jo 


00 



1 



AXg^'Kcos ncp' sin stp-sin nfip' cos wpjdcp' 
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Ay x (qp) = b Q + i ^ ± ^ Aii( < P’) Bin n(cp* - qp) dcp* - 


1 V * N 2 P2jt 

« / "n n / Ax 2^') sin n(qp’- qp) dqp* ( 

1 * • * ^0 

Nov let f (qp) be a function that can be developed in a Fourier 
series for 0 = cp= 2 a. Then 

J' f(qp') cot d( P* ~ y " ^ f(qp') sin n(qp'-qp) dqp' = 

® m=l 


0 

(A7) 


Hence, 




Ax 1 (qp' ) sin n(qp 


• - qp)J - 


bi^cp*) |j| cot .C^g . . r ?p) _ sin n(qp' -qp)jj dc P' - 


op gj( 

* 22 ^ : q -n Sg “2 (q>,) 8l " n(<p ' • ^ 


dqp' 


or 


Ayx^) =^ 0 -^ u Ax 1 (cp')|| cot +^i- sin n(cp'-cp)J dcp’ 


\ ' f AXgCqp') sin n(qp'-cp) dcp* 


(A8) 


Inasmuch as the series of equation (A8) are uniformly convergent as 
are the series of equations (5). the summation and integration may 
be interchanged in equation (A8) to yield equation (5). 
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APPENDIX 8 


THE NUMERICAL EVALUATION OF THE CARTESIAN MAPPING FUNCTION 

The det erminat ion of the functions Ax^ and AX£ from the 
given functions Ay^ and Ay 2 was based in this paper on numerical 
Integration of the first two equations (5) . The equations for Ax-^ 
and AXg; when the constant a q has been set equal to zero, are 


axjfo) - i 


r2jt 


Ay^cp') 


5 cot ^ + \ 2 - S - -- Bin n(<p'-<P) 


V ^ 


dcp'- 


> (Bl) 


dcp'+ 


Ay 2 (cp') ~^ 2 n sin n ( c * >, - ( P) 4C P' 

Ax 2 (cp)=-~ f Ay (cp* ) | cot + \ sin n(cp'-cp) 

JO 2 c Z — f 1-q 

L 1 

p2jt o° n 

j-j Ay^cp') ^ 2 %n sin n ( <p,_<T> ) 4( P' 

Jo 1 


If the range of cp is divided into 2n equal Intervals whose length 
is 5, if the values of Ay are given at the end points of the 
intervals, and if Ax is desired at the same points, approximate 
integration will yield expressions of the following form: 

2n-l 2n-l 

Ail(<P) = 1 (cjc + dfc) Ay x (cp + MB) + > ^ Ay 2 (<P + MB) 


M=0 


k=0 


2n-l 2n-l 

Ax 2 (cp) ■- >' ^ (c k + d k ) Ay 2 (cp + k6) - *> ' ej^. Ay^cp + MB) 


M=0 


k=0 


(B2) 


The values of Cj^ have been calculated in reference 1 by means of 
Simpson's rule and other simplifications for use with the CMF of 
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simply connected regions. The values of d^ and e^ may he sim- 
ilarly obtained. The value of c k as calculated in reference 1 is 


c 0 = 0 


6 . 5 

°1 - S cot 2 


'N 


8 + sin 8 
2a sin 8 


-8 _ A 8 8 + sin 8 [ 

c 2n-l = 6rt Cot 2 “ 2it sin 6 r 

c k s ^ cot T ^ 


(B3) 


= H COt T (k eV6n) J 


In the present paper, because the number of intervals was an inte- 
gral multiple of 6, Weddle's rule was used for the evaluation of d k 
a nd eif. 

The values of c^ are given in table XV for the cases of 
2n = 24 2n = 48. The values of djj. and e^ contain the 

parameter q. Hence, these coefficients must be evaluated anew for 
each approximation. 
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TABLE I - OEDINATES OF 12-PERCENT THICK AIRFOIL 
[From table 1 of reference 1] 


Station 

(percent 

chord 

from 

nose) 

Ordinate 


Station 

(percent 

chord 

from 

nose) 

Ordinate 

0 

0 


50 

5.880 

1.25 

1.425 


55 

5.540 

2.5 

1.900 


60 

5.025 

5 

2.585 


65 

4.415 

10 

3.540 


70 

3.750 

15 

4.250 


75 

3.060 

20 

4.820 


80 

2.350 

25 

5.295 


85 

1.685 

30 

5.655 


90 

1.060 

35 

5.900 


95 

.510 

40 

6.000 


97.5 

.260 

45 

6.010 


100 

0 




TABLE II - VELOCITY DISTRIBUTION AND CARTESIAN MAPPING FUNCTION FOR AIRFOIL AT ANGLE OF ATTACK OF 0° 
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TABLE III - ■VELOCITY DISTRIBUTION AND CARTESIAN MAPPING 
FUNCTION FOR AIRFOIL AT ANGLE OF ATTACK OF 4° 

[q. = 0.2041; q> N = 94°; <P T = 274°; q > 2 = 89.91°; T = -0.0063] 



Airfoil 

Channel walls 

(deg) 

4l i 

v cl 

r(cp) 

Az l 

&AX-l 

dAy x 

4x 2 

v c2 

< 

Ax 2 

dAxg 

dAy 2 

V 

dqp 

dcp 

y 

dcp 

dqp 

0 X 7.5 

0.0152 

1.0735 

-0.02940 

0.01005 

0.0376 

0.0033 

0.0236 

0.9730 

0 

0.00591 

0.0114 

0 

1 

-.1010 

1.0709 

-.02840 

.01511 

.0372 

.0110 

-.1917 

.9719 

0 

.00739 

.0109 

0 

2 

-.2168 

1.0617 

-.02645 

.01983 

.0355 

.0185 

-.3850 

.9703 

0 

.00875 

.0101 

0 

3 

-.3308 

1.0419 

-.02380 

.02415 

.0316 

.0244 

-.5838 

.9686 

0 

.01000 

.0090 

0 

4 

-.4417 

1.0133 

-.02028 

.02795 

.0260 

.0291 

-.7920 

.9670 

0 

.01107 

.0070 

0 

5 

-.5486 

.9822 

-.01622 

.03095 

.0203 

.0330 

-1.0134 

.9659 

0 

.01198 

.0062 

0 

6 

-.6493 

.9464 

-.01182 

.03317 

.0145 

.0356 

-1.2549 

.9657 

0 

.01267 

.0045 

0 

7 

-.7415 

.9036 

-.00720 

.03465 

.0088 

.0369 

-1.5261 

.9668 

0 

.01316 

.0027 

0 

8 

-.8234 

.8247 

-.00242 

.03528 

-.0013 

.0355 

-1.8436 

.9697 

0 

.01340 

.0010 

0 

9 

-.8928 

.7729 

.00212 

.03502 

-.0020 

.0313 

-2.2367 

.9743 

0 

.01342 

-.0009 

0 

10 

-.9448 

.6579 

.00598 

.03464 

-.0030 

.0303 

-2.7752 

.9810 

0 

.01319 

-.0027 

0 

11 

-.9776 

.3130 

.01035 

.03418 

-.0082 

.0347 

-3.6780 

.9896 

0 

.01272 

-.0045 

0 

12 

-.9950 

.5485 

.01525 

.03241 

-.0194 

.0338 

-9.2779 

1.0001 

0 

.01202 

-.0062 

0 

13 

-.9963 

1.7756 

.01925 

.02947 

-.0229 

.0299 

-3.6554 

1.0121 

0 

.01109 

-.0079 

0 

14 

-.9777 

1.6934 

.02312 

.02634 

-.0252 

.0250 

-2.7736 

1.0249 

0 

.00997 

-.0094 

0 

15 

-.9407 

1.5167 

.02600 

.02294 

-.0275 

.0210 

-2.2462 

1.0379 

0 

.00864 

-.0106 

0 

16 

-.8875 

1.4150 

.02862 

.01915 

-.0295 

.0183 

-1.8613 

1.0506 

0 

.00718 

-.0118 

0 

17 

-.8188 

1,3589 

.03075 

.01518 

-.0315 

.0155 

-1.5513 

1.0624 

0 

.00555 

-.0127 

0 

18 

-.7375 

1.3249 

.03252 

.01096 

-.0333 

.0121 

-1.2863 

1.0731 

0 

.00384 

-.0135 

0 

19 

-.6456 

1.3042 

.03390 

.00652 

-.0350 

.0085 

-1.0503 

1.0817 

0 

.00203 

-.0140 

0 

20 

-.5455 

1.2910 

.03475 

.00184 

-.0365 

.0048 

-.8330 

1.0883 

0 

.00021 

-.0142 

0 

21 

-.4391 

1.2843 

.03522 

-.00305 

-.0380 

.0005 

-.6288 

1.0926 

0 

-.00164 

-.0141 

0 

22 

-.3287 

1.2803 

.03502 

-.00824 

-.0392 

-.0047 

-.4326 

1.0948 

0 

-.00343 

-.0138 

0 

23 

-.2152 

1.2790 

.03402 

-.01349 

-.0402 

-.0101 

-.2414 

1.0944 

0 

-.00518 

-.0131 

0 

24 

-.1004 

1.2772 

.03235 

-.01883 

-.0406 

-.0162 

-.0272 

1.0918 

0 

-.00680 

-.0120 

0 

25 

.0143 

1.2652 

.02970 

-.02423 

-.0390 

-.0250 

.1375 

1.0874 

0 

-.00831 

-.0107 

0 

26 

.1287 

1.2392 

.02558 

-.02926 

-.0350 

-.0364 

.3303 

1.0814 

0 

-.00963 

-.0092 

0 

27 

.2439 

1.1935 

.02038 

-.03325 

-.0263 

-.04281 

.5290 

1.0745 

0 

-.01079 

-.0077 

0 

28 

.3592 

1.1455 

.01438 

-.03590 

-.0160 

-.0456 

.7370 

1.0667 

0 

-.01171 

-.0061 

0 

29 

.4732 

1.0927 

.00828 

-.03712 

-.0041 

-.0462 

.9585 

1.0576 

0 

-.01243 

-.0042 

0 

30 

.5834 

1.0475 

.00215 

-.03695 

.0065 

-.0441 

1.2000 

1.0491 

0 

-.01290 

-.0027 

0 

31 

.6871 

1.0093 

-.00325 

- .03558 

.0153 

-.0395 

1.4711 

1.0400 

0 

-.01314 

-.0010; 

0 

32 

.7812 

.9805 

-.00815 

-.03317 

.0216 

-.0314 

1.7876 

1.0308 

0 

-.01315 

.0008 

0 

33 

.8621 

.9564 

-.01178 

-.03007 

.0253 

-.0240 

2.1793 

1.0221 

0 

-.01294 

.0023 

0 

34 

.9260 

.9328 

-.01450 

-.02674 

.0269 

-.0175 

2.7132 

1.0139 

0 

-.01252 

.0039 

0 

35 

.9708 

.9223 

-.01642 

-.02332 

.0254 

-.0110 

3.6019 

1.0066 

0 

-.01190 

.0053 

0 

36 

.9948 

.9246 

-.01742 

-.01996 

.0225 

-.0043 

9.2833 

1.0001 

0 

-.01112 

.0066 

0 

37 

.9953 

.8741 

-.01750 

-.01723 

.0185 

0 

3.6385 

.9946 

0 

-.01016 

.0078 

0 

38 

.9723 

.8882 

-.01760 

-.01526 

.0140 

-.0012 

2.7416 

.9900 

0 

-.00907 

.0088 

0 

39 

.9280 

.9050 

-.01792 

-.01370 

.0115 

-.0030 

2.2091 

.9864 

0 

-.00785 

.0096 

0 

40 

.8653 

.9190 

-.01850 

-.01229 

.0102 

-.0051 

1.8208 

.9836 

0 

-.00654 

.0104 

0 

41 

.7861 

.9307 

-.01932 

-.01099 

.0099 

-.0080 

1.5084 

.9814 

0 

-.00511 

.0110 

0 

42 

.6932 

.9427 

- .02050 

-.00968 

.0107 

-.0116 

1.2411 

.9800 

0 

-.00362 

.0116 

0 

43 

.5898 

.9608 

-.02218 

-.00813 

.0136 

-.0141 

1.0030 

.9789 

0 

-.00206 

.0120 

0 

44 

.4796 

.9876 

-.02425 

-.00594 

.0190 

-.0148 

.7842 

.9780 

0 

-.00048 

.0123 

0 

45 

.3653 

1.0142 

-.02620 

-.00300 

.0248 

-.0146 

.5784 

.9772 

0 

.00116 

.0124 

0 

46 

.2485 

1.0414 

-.02788 

.00061 

.0308 

-.0127 

.3809 

.9762 

0 

.00276 

.0123 

0 

47 

.1314 

1.0669 

-.02920 

.00501 

.0362 

-.0064 

.1887 

.9748 

0 

.00437 

.0120 

0 
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TABLE IV - COEFFICIENTS FOR CALCULATION OF CARTESIAN 
MAPPING FUNCTION FOR SINGLE CONTOUR 
(a) 24-point scheme 


k 

c k 

k 

c k 

0 

0 

12 

0 

1 

.42564 

13 

-.00366 

2 

.20734 

14 

-.01489 

3 

.06706 

15 

-.01151 

4 

.09623 

16 

-.03208 

5 

.03620 

17 

-.02131 

6 

.05556 

18 

-.05556 

7 

.02131 

19 

-.03620 

8 

.03208 

20 

-.09623 

9 

.01151 

21 

-.06706 

10 

.01489 

22 

-.20734 

11 

.00366 

23 

-.42564 

(b) 48-point 

scheme 

k 

c k 

k 

c k 

0 

0 

24 

0 

1 

.42470 

25 

-.00091 

2 

.21099 

26 

-.00366 

3 

.06982 

27 

-.00276 

4 

.10367 

28 

-.00744 

5 

.04092 

29 

-.00472 

6 

.06706 

30 

-.01151 

7 

.02816 

31 

-.00685 

8 

.04811 

32 

-.01604 

9 

.02079 

33 

-.00928 

10 

.03620 

34 

-.02132 

11 

.01584 

35 

-.01218 

12 

.02778 

36 

-.02778 

13 

.01218 

37 

-.01584 

14 

.02132 

38 

-.03620 

15 

.00928 

39 

-.02079 

16 

.01604 

40 

-.04811 

17 

.00685 

41 

-.02816 

18 

.01151 

42 

-.06706 

19 

.00472 

43 

-.04092 

20 

.00744 

44 

-.10367 

21 

.00276 

45 

-.06982 

22 

.00366 

46 

-.21099 

23 

.00091 

47 

-.42470 
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Channel wall H * — >• 



Center line of channel 


«♦— H Channel wall 


4-plane 



Figure 1.- Transformation of flat plate and 
channel into two concentric circles. 



Channel wall in z-plane 
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Figure 3, - Velocity distribution on airfoil in channel at angle of attack of 0° and velocity distribution on c hann el walls. 





Velocity distribution, V c ]/ v 
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Figure 4. - Comparison of velocity distributions in channel obtained by two methods . 



wall 
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Figure 5, - Velocity distribution on airfoil in channel at angle of attack of 4° and velocity distribution on channel vails 
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Figure 6. - Comparison of velocities on airfoil in free stream and on airfoil in channel 

for a = 4°. 







Local velocity correction 
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("b) Lower surface. 

Figure 8. - Velocity corrections for 12-percent-thick airfoil, a = 4°, c/h = 0.5. 
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